Below, we first show an aggregation of the latest version of FNGS for our respective post-processing pipelines:

  require(Discriminability)
## Loading required package: Discriminability
  signal_overall <- c()
  ids_overall <- c()
  cpac_graphs_overall <- c()
  cpac_ids_overall <- c()

FNGS

BNU1

  inpath <- 'C:/Users/ebrid/Documents/R/fngs_merge/BNU1/roi_timeseries/desikan-2mm/'
  tsnames <- list.files(inpath, pattern="\\.rds", full.names=TRUE)
  scan_pos = 3
  tsobj <- open_timeseries(tsnames, sub_pos=scan_pos, exclude = TRUE)
## [1] "opening timeseries..."
  signal <- tsobj[[1]]
  ids <- tsobj[[3]]
  
  signal_overall <- c(signal_overall, signal)
  ids_overall <- c(ids_overall, tsobj$dataset)

  cpac_inpath <- 'C:/Users/ebrid/Documents/R/CPAC_results/BNU1/FSL_nff_nsc_ngs_des/'
  cpac_tsnames <- list.files(cpac_inpath, pattern="\\.graphml", full.names=TRUE)
  scan_pos = 2
  cpac_graobj <- open_graphs(cpac_tsnames, sub_pos=scan_pos)
## [1] "opening graphs..."
## [1] 1
## [1] 2
## [1] 3
## [1] 4
## [1] 5
## [1] 6
## [1] 7
## [1] 8
## [1] 9
## [1] 10
## [1] 11
## [1] 12
## [1] 13
## [1] 14
## [1] 15
## [1] 16
## [1] 17
## [1] 18
## [1] 19
## [1] 20
## [1] 21
## [1] 22
## [1] 23
## [1] 24
## [1] 25
## [1] 26
## [1] 27
## [1] 28
## [1] 29
## [1] 30
## [1] 31
## [1] 32
## [1] 33
## [1] 34
## [1] 35
## [1] 36
## [1] 37
## [1] 38
## [1] 39
## [1] 40
## [1] 41
## [1] 42
## [1] 43
## [1] 44
## [1] 45
## [1] 46
## [1] 47
## [1] 48
## [1] 49
## [1] 50
## [1] 51
## [1] 52
## [1] 53
## [1] 54
## [1] 55
## [1] 56
## [1] 57
## [1] 58
## [1] 59
## [1] 60
## [1] 61
## [1] 62
## [1] 63
## [1] 64
## [1] 65
## [1] 66
## [1] 67
## [1] 68
## [1] 69
## [1] 70
## [1] 71
## [1] 72
## [1] 73
## [1] 74
## [1] 75
## [1] 76
## [1] 77
## [1] 78
## [1] 79
## [1] 80
## [1] 81
## [1] 82
## [1] 83
## [1] 84
## [1] 85
## [1] 86
## [1] 87
## [1] 88
## [1] 89
## [1] 90
## [1] 91
## [1] 92
## [1] 93
## [1] 94
## [1] 95
## [1] 96
## [1] 97
## [1] 98
## [1] 99
## [1] 100
  cpac_graphs <- cpac_graobj[[1]]
  cpac_ids <- cpac_graobj[[3]]
  
  cpac_graphs_overall <- c(cpac_graphs_overall, cpac_graphs)
  ids_overall <- c(ids_overall, cpac_graobj$dataset)

Sanity Check

Timeseries

  ts1 = signal[[1]]
  title = paste('Timeseries plot for subject', ids[1])
  print(plot_timeseries(ts1, title=title))
## Warning: package 'reshape2' was built under R version 3.3.3
## Warning: package 'ggplot2' was built under R version 3.3.3

#### FNGS Correlation Matrix

  title = paste('Correlation plot for subject', ids[1])
  print(plot_square(abs(cor(ts1)), title=title, legend="corr"))

CPAC Correlation Matrix

  title = paste('Correlation plot for subject', ids[1])
  print(plot_square(cpac_graphs[[1]], title=title, legend="corr"))

FNGS Raw Correaltion Graphs

  td <- time_discr(signal, ids, rank=FALSE)
## [1] "computing hellinger distances and densities..."

  print(td$combinedplot)
## NULL

FNGS Ranked Correlation Graphs

  tdr <- time_discr(signal, ids, rank=TRUE)
## [1] "computing hellinger distances and densities..."

  print(tdr$combinedplot)
## NULL

CPAC Raw Correlation Graphs

  cpac_td <- time_discr(cpac_graphs, cpac_ids, rank=FALSE, graphs = TRUE)
## [1] "computing hellinger distances and densities..."

  print(cpac_td$combinedplot)
## NULL

CPAC Ranked Correlation Graphs

  cpac_tdr <- time_discr(cpac_graphs, cpac_ids, rank=TRUE, graphs = TRUE)
## [1] "computing hellinger distances and densities..."

  print(cpac_tdr$combinedplot)
## NULL
  collection <- data.frame(dataset='BNU1',
                                             postprocessing=c('fngs_corr', 'fngs_rcorr',
                                                              'cpac_corr', 'cpac_rcorr'),
                           discr=c(td$d, tdr$d, cpac_td$d, cpac_tdr$d))

DC1

  inpath <- 'C:/Users/ebrid/Documents/R/fngs_merge/DC1/roi_timeseries/desikan-2mm/'
  tsnames <- list.files(inpath, pattern="\\.rds", full.names=TRUE)
  scan_pos = 3
  tsobj <- open_timeseries(tsnames, sub_pos=scan_pos, exclude = TRUE)
## [1] "opening timeseries..."
  signal <- tsobj[[1]]
  ids <- tsobj[[3]]
  
  signal_overall <- c(signal_overall, signal)
  ids_overall <- c(ids_overall, tsobj$dataset)

  cpac_inpath <- 'C:/Users/ebrid/Documents/R/CPAC_results/DC1/FSL_nff_nsc_ngs_des/'
  cpac_tsnames <- list.files(cpac_inpath, pattern="\\.graphml", full.names=TRUE)
  scan_pos = 2
  cpac_graobj <- open_graphs(cpac_tsnames, sub_pos=scan_pos)
## [1] "opening graphs..."
## [1] 1
## [1] 2
## [1] 3
## [1] 4
## [1] 5
## [1] 6
## [1] 7
## [1] 8
## [1] 9
## [1] 10
## [1] 11
## [1] 12
## [1] 13
## [1] 14
## [1] 15
## [1] 16
## [1] 17
## [1] 18
## [1] 19
## [1] 20
## [1] 21
## [1] 22
## [1] 23
## [1] 24
## [1] 25
## [1] 26
## [1] 27
## [1] 28
## [1] 29
## [1] 30
## [1] 31
## [1] 32
## [1] 33
## [1] 34
## [1] 35
## [1] 36
## [1] 37
## [1] 38
## [1] 39
## [1] 40
## [1] 41
## [1] 42
## [1] 43
## [1] 44
## [1] 45
## [1] 46
## [1] 47
## [1] 48
## [1] 49
## [1] 50
## [1] 51
## [1] 52
## [1] 53
## [1] 54
## [1] 55
## [1] 56
## [1] 57
## [1] 58
## [1] 59
## [1] 60
## [1] 61
## [1] 62
## [1] 63
## [1] 64
## [1] 65
## [1] 66
## [1] 67
## [1] 68
## [1] 69
## [1] 70
## [1] 71
## [1] 72
## [1] 73
## [1] 74
## [1] 75
## [1] 76
## [1] 77
## [1] 78
## [1] 79
## [1] 80
## [1] 81
## [1] 82
## [1] 83
## [1] 84
## [1] 85
## [1] 86
## [1] 87
## [1] 88
## [1] 89
## [1] 90
## [1] 91
## [1] 92
## [1] 93
## [1] 94
## [1] 95
## [1] 96
## [1] 97
## [1] 98
## [1] 99
## [1] 100
## [1] 101
## [1] 102
## [1] 103
## [1] 104
## [1] 105
## [1] 106
## [1] 107
## [1] 108
## [1] 109
## [1] 110
## [1] 111
## [1] 112
## [1] 113
## [1] 114
## [1] 115
## [1] 116
## [1] 117
## [1] 118
## [1] 119
## [1] 120
## [1] 121
## [1] 122
## [1] 123
## [1] 124
## [1] 125
## [1] 126
## [1] 127
## [1] 128
## [1] 129
## [1] 130
## [1] 131
## [1] 132
## [1] 133
## [1] 134
## [1] 135
## [1] 136
## [1] 137
## [1] 138
## [1] 139
## [1] 140
## [1] 141
## [1] 142
## [1] 143
## [1] 144
## [1] 145
## [1] 146
## [1] 147
## [1] 148
## [1] 149
## [1] 150
## [1] 151
## [1] 152
## [1] 153
## [1] 154
## [1] 155
## [1] 156
## [1] 157
## [1] 158
## [1] 159
## [1] 160
## [1] 161
## [1] 162
## [1] 163
## [1] 164
## [1] 165
## [1] 166
## [1] 167
## [1] 168
## [1] 169
## [1] 170
## [1] 171
## [1] 172
## [1] 173
## [1] 174
## [1] 175
## [1] 176
## [1] 177
## [1] 178
## [1] 179
## [1] 180
## [1] 181
## [1] 182
## [1] 183
## [1] 184
## [1] 185
## [1] 186
## [1] 187
## [1] 188
## [1] 189
## [1] 190
## [1] 191
## [1] 192
## [1] 193
## [1] 194
## [1] 195
## [1] 196
## [1] 197
## [1] 198
## [1] 199
## [1] 200
## [1] 201
## [1] 202
## [1] 203
## [1] 204
## [1] 205
## [1] 206
## [1] 207
## [1] 208
## [1] 209
## [1] 210
## [1] 211
## [1] 212
## [1] 213
## [1] 214
## [1] 215
## [1] 216
## [1] 217
## [1] 218
## [1] 219
## [1] 220
## [1] 221
## [1] 222
## [1] 223
## [1] 224
## [1] 225
## [1] 226
## [1] 227
## [1] 228
## [1] 229
## [1] 230
## [1] 231
## [1] 232
## [1] 233
## [1] 234
## [1] 235
## [1] 236
## [1] 237
## [1] 238
## [1] 239
## [1] 240
## [1] 241
## [1] 242
## [1] 243
## [1] 244
  cpac_graphs <- cpac_graobj[[1]]
  cpac_ids <- cpac_graobj[[3]]
  
  cpac_graphs_overall <- c(cpac_graphs_overall, cpac_graphs)
  ids_overall <- c(ids_overall, cpac_graobj$dataset)

Sanity Check

Timeseries

  ts1 = signal[[1]]
  title = paste('Timeseries plot for subject', ids[1])
  print(plot_timeseries(ts1, title=title))

#### FNGS Correlation Matrix

  title = paste('Correlation plot for subject', ids[1])
  print(plot_square(abs(cor(ts1)), title=title, legend="corr"))

CPAC Correlation Matrix

  title = paste('Correlation plot for subject', ids[1])
  print(plot_square(cpac_graphs[[1]], title=title, legend="corr"))

FNGS Raw Correaltion Graphs

  td <- time_discr(signal, ids, rank=FALSE)
## [1] "computing hellinger distances and densities..."

  print(td$combinedplot)
## NULL

FNGS Ranked Correlation Graphs

  tdr <- time_discr(signal, ids, rank=TRUE)
## [1] "computing hellinger distances and densities..."

  print(tdr$combinedplot)
## NULL

CPAC Raw Correlation Graphs

  cpac_td <- time_discr(cpac_graphs, cpac_ids, rank=FALSE, graphs = TRUE)
## [1] "computing hellinger distances and densities..."

  print(cpac_td$combinedplot)
## NULL

CPAC Ranked Correlation Graphs

  cpac_tdr <- time_discr(cpac_graphs, cpac_ids, rank=TRUE, graphs = TRUE)
## [1] "computing hellinger distances and densities..."

  print(cpac_tdr$combinedplot)
## NULL
  collection <- rbind(collection, data.frame(dataset='DC1',
                                             postprocessing=c('fngs_corr', 'fngs_rcorr',
                                                              'cpac_corr', 'cpac_rcorr'),
                           discr=c(td$d, tdr$d, cpac_td$d, cpac_tdr$d)))

HNU1

  inpath <- 'C:/Users/ebrid/Documents/R/fngs_merge/HNU1/roi_timeseries/desikan-2mm/'
  tsnames <- list.files(inpath, pattern="\\.rds", full.names=TRUE)
  scan_pos = 3
  tsobj <- open_timeseries(tsnames, sub_pos=scan_pos, exclude = TRUE)
  signal <- tsobj[[1]]
  ids <- tsobj[[3]]
  
  signal_overall <- c(signal_overall, signal)
  ids_overall <- c(ids_overall, tsobj$dataset)

  cpac_inpath <- 'C:/Users/ebrid/Documents/R/CPAC_results/HNU1/FSL_nff_nsc_ngs_des/'
  cpac_tsnames <- list.files(cpac_inpath, pattern="\\.graphml", full.names=TRUE)
  scan_pos = 2
  cpac_graobj <- open_graphs(cpac_tsnames, sub_pos=scan_pos)
  cpac_graphs <- cpac_graobj[[1]]
  cpac_ids <- cpac_graobj[[3]]
  
  cpac_graphs_overall <- c(cpac_graphs_overall, cpac_graphs)
  ids_overall <- c(ids_overall, cpac_graobj$dataset)

Sanity Check

Timeseries

  ts1 = signal[[1]]
  title = paste('Timeseries plot for subject', ids[1])
  print(plot_timeseries(ts1, title=title))

FNGS Correlation Matrix

  title = paste('Correlation plot for subject', ids[1])
  print(plot_square(abs(cor(ts1)), title=title, legend="corr"))

CPAC Correlation Matrix

  title = paste('Correlation plot for subject', ids[1])
  print(plot_square(cpac_graphs[[1]], title=title, legend="corr"))

FNGS Raw Correaltion Graphs

  td <- time_discr(signal, ids, rank=FALSE)
  print(td$combinedplot)

FNGS Ranked Correlation Graphs

  tdr <- time_discr(signal, ids, rank=TRUE)
  print(tdr$combinedplot)

CPAC Raw Correlation Graphs

  cpac_td <- time_discr(cpac_graphs, cpac_ids, rank=FALSE, graphs = TRUE)
  print(cpac_td$combinedplot)

CPAC Ranked Correlation Graphs

  cpac_tdr <- time_discr(cpac_graphs, cpac_ids, rank=TRUE, graphs = TRUE)
  print(cpac_tdr$combinedplot)
  collection <- rbind(collection, data.frame(dataset='HNU1',
                                             postprocessing=c('fngs_corr', 'fngs_rcorr',
                                                              'cpac_corr', 'cpac_rcorr'),
                           discr=c(td$d, tdr$d, cpac_td$d, cpac_tdr$d)))

SWU4

  inpath <- 'C:/Users/ebrid/Documents/R/fngs_merge/SWU4/roi_timeseries/desikan-2mm/'
  tsnames <- list.files(inpath, pattern="\\.rds", full.names=TRUE)
  scan_pos = 3
  tsobj <- open_timeseries(tsnames, sub_pos=scan_pos, exclude = TRUE)
  signal <- tsobj[[1]]
  ids <- tsobj[[3]]
  
  signal_overall <- c(signal_overall, signal)
  ids_overall <- c(ids_overall, tsobj$dataset)

  cpac_inpath <- 'C:/Users/ebrid/Documents/R/CPAC_results/SWU4/FSL_nff_nsc_ngs_des/'
  cpac_tsnames <- list.files(cpac_inpath, pattern="\\.graphml", full.names=TRUE)
  scan_pos = 2
  cpac_graobj <- open_graphs(cpac_tsnames, sub_pos=scan_pos)
  cpac_graphs <- cpac_graobj[[1]]
  cpac_ids <- cpac_graobj[[3]]
  
  cpac_graphs_overall <- c(cpac_graphs_overall, cpac_graphs)
  ids_overall <- c(ids_overall, cpac_graobj$dataset)

Sanity Check

Timeseries

  ts1 = signal[[1]]
  title = paste('Timeseries plot for subject', ids[1])
  print(plot_timeseries(ts1, title=title))

FNGS Correlation Matrix

  title = paste('Correlation plot for subject', ids[1])
  print(plot_square(abs(cor(ts1)), title=title, legend="corr"))

CPAC Correlation Matrix

  title = paste('Correlation plot for subject', ids[1])
  print(plot_square(cpac_graphs[[1]], title=title, legend="corr"))

FNGS Raw Correaltion Graphs

  td <- time_discr(signal, ids, rank=FALSE)
  print(td$combinedplot)

FNGS Ranked Correlation Graphs

  tdr <- time_discr(signal, ids, rank=TRUE)
  print(tdr$combinedplot)

CPAC Raw Correlation Graphs

  cpac_td <- time_discr(cpac_graphs, cpac_ids, rank=FALSE, graphs = TRUE)
  print(cpac_td$combinedplot)

CPAC Ranked Correlation Graphs

  cpac_tdr <- time_discr(cpac_graphs, cpac_ids, rank=TRUE, graphs = TRUE)
  print(cpac_tdr$combinedplot)

SWU4

  inpath <- 'C:/Users/ebrid/Documents/R/fngs_merge/BNU2/roi_timeseries/desikan-2mm/'
  tsnames <- list.files(inpath, pattern="\\.rds", full.names=TRUE)
  scan_pos = 3
  tsobj <- open_timeseries(tsnames, sub_pos=scan_pos, exclude = TRUE)
  signal <- tsobj[[1]]
  ids <- tsobj[[3]]
  
  signal_overall <- c(signal_overall, signal)
  ids_overall <- c(ids_overall, tsobj$dataset)

  cpac_inpath <- 'C:/Users/ebrid/Documents/R/CPAC_results/BNU2/FSL_nff_nsc_ngs_des/'
  cpac_tsnames <- list.files(cpac_inpath, pattern="\\.graphml", full.names=TRUE)
  scan_pos = 2
  cpac_graobj <- open_graphs(cpac_tsnames, sub_pos=scan_pos)
  cpac_graphs <- cpac_graobj[[1]]
  cpac_ids <- cpac_graobj[[3]]
  
  cpac_graphs_overall <- c(cpac_graphs_overall, cpac_graphs)
  ids_overall <- c(ids_overall, cpac_graobj$dataset)

Sanity Check

Timeseries

  ts1 = signal[[1]]
  title = paste('Timeseries plot for subject', ids[1])
  print(plot_timeseries(ts1, title=title))

FNGS Correlation Matrix

  title = paste('Correlation plot for subject', ids[1])
  print(plot_square(abs(cor(ts1)), title=title, legend="corr"))

CPAC Correlation Matrix

  title = paste('Correlation plot for subject', ids[1])
  print(plot_square(cpac_graphs[[1]], title=title, legend="corr"))

FNGS Raw Correaltion Graphs

  td <- time_discr(signal, ids, rank=FALSE)
  print(td$combinedplot)

FNGS Ranked Correlation Graphs

  tdr <- time_discr(signal, ids, rank=TRUE)
  print(tdr$combinedplot)

CPAC Raw Correlation Graphs

  cpac_td <- time_discr(cpac_graphs, cpac_ids, rank=FALSE, graphs = TRUE)
  print(cpac_td$combinedplot)

CPAC Ranked Correlation Graphs

  cpac_tdr <- time_discr(cpac_graphs, cpac_ids, rank=TRUE, graphs = TRUE)
  print(cpac_tdr$combinedplot)
  collection <- rbind(collection, data.frame(dataset='HNU1',
                                             postprocessing=c('fngs_corr', 'fngs_rcorr',
                                                              'cpac_corr', 'cpac_rcorr'),
                           discr=c(td$d, tdr$d, cpac_td$d, cpac_tdr$d)))
  collection$dataset <- factor(collection$dataset)
  collection$postprocessing <- factor(collection$postprocessing)
  summaryfig_fng <- ggplot(collection, aes(x=postprocessing, y=discr, color=dataset)) +
    geom_point(size=4) +
    ggtitle('Comparing Discriminability over 4 reference datasets using different postprocessing techniques')
  print(summaryfig_fng)

  summaryfig_fng <- ggplot(collection, aes(x=postprocessing, y=discr, fill=postprocessing)) +
    geom_violin() +
    geom_boxplot(width=0.1) +
    ggtitle('Comparing Discriminability over 4 reference datasets using different postprocessing techniques')
  print(summaryfig_fng)